Ground-state structure of the hydrogen double vacancy on Pd(lll) 
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We determine the ground-state structure of a double vacancy in a hydrogen monolayer on the 
Pd(lll) surface. We represent the double vacancy as a triple vacancy containing one additional 
hydrogen atom. The potential-energy surface for a hydrogen atom moving in the triple vacancy is 
obtained by density-functional theory, and the wave function of the fully quantum hydrogen atom 
is obtained by solving the Schrodinger equation. We find that an H atom in a divacancy defect 
experiences significant quantum effects, and that the ground-state wave function is centered at the 
hep site rather than the fee site normally occupied by H atoms on Pd(lll). Our results agree well 
with scanning tunneling microscopy images. 

PACS numbers: 31.15.es, 61.72.jd, 68.43.Fg, 



I. INTRODUCTION 

Understanding the fundamental mechanisms that gov- 
ern the behavior of hydrogen on metal surfaces is im- 
portant to new technologies involving hydrogen produc- 
tion, storage, and energy conversioni^^i^i^ . Because of 
its light mass, hydrogen can manifest uniquely quantum 
effects not seen for other elements. These effects can 
have striking consequences when the dimensionality of 
the allowed motion is constrained by adsorption on a 
surface^' -iii^i^ii^. When a metal surface such as palla- 
dium is covered by a nearly complete monolayer of hy- 
drogen, small clusters of vacancies in the hydrogen layer 
serve as active sites for the dissociative adsorption of ad- 
ditional H2 moleculesiiii^. Despite their central role in 
catalysis, the exact nature of these vacancy defects is not 
well understood. 

It has been generally accepted that the dissociative 
adsorption of the diatomic molecule H2 is to require at 
least two adjacent and empty atomic adsorption sites 
(or vacancies) based on early studies of Langmuir^^^^^. 
Recently, intriguing observations about hydrogen diva- 
cancies on palladium were made using scanning tunnel- 
ing microscopy (STM). Mitsui and co-worker^i^^i^ re- 
ported that hydrogen molecules impinging on an al- 
most H-saturated Pd(lll) surface did not adsorb in 
two- vacancy sites. Rather, aggregates of three or more 
hydrogen vacancies were required for efficient dissocia- 
tion of H2 molecules. These findings led to speculations 
that the standard description from Langmuir adsorption 
kineticsi^ might be too simple to explain the dissociative 
adsorption of hydrogen on Pd(lll)i2iii^ii^ii^. 

However, GroB and Dianat have recently shown, us- 
ing ah initio molecular dynamics (AIMD) simulations, 
that no change in the original Langmuir picture is 
warranted^^. Based on more than 4000 AIMD trajecto- 
ries, GroB et at. unambiguously demonstrated that a di- 
vacancy is still sufficient to dissociate hydrogen molecules 



provided the kinetic energy of the incident molecules 
is large enough to overcome the relatively small en- 
ergy barrier. For an initial kinetic energy of 0.02 eV, 
roughly corresponding to the kinetic energy of the II2 
molecules in Mitsui's experiment, no dissociative ad- 
sorption events were observed on the almost hydrogen- 
covered Pd(lll) surface, thus confirming the experimen- 
tal findings. When the initial kinetic energy was in- 
creased to 0.1 eV so that the small barrier found for 
H2 on Pd(lll) at high coverage^^ can be overcome, a 
nonvanishing adsorption probability was observed. 

To further improve our understanding of dissociative 
adsorption of H2 on Pd(lll), an accurate description of 
the ground-state structure of the double vacancy in a hy- 
drogen monolayer on the surface is desirable. Mitsui et 
^^14,15 found that divacancies have a triangular STM im- 
age extending over three neighboring sites, instead of a 
linear image over two neighboring sites. In this paper, by 
determining the ground-state structure of the divacancy 
in a hydrogen monolayer on Pd(lll), we show theoret- 
ically that the quantum nature of the H atom accounts 
for the observed STM imagery. Our findings indicate 
that the quantum wave nature of hydrogen motion pro- 
vides a simple and elegant foundation to the traditional 
classical explanations for such observationsi^ii^ii^, and 
suggest that the quantum nature of hydrogen may play 
a surprisingly prominent role in other similarly confined 
systems. 



II. METHODS 

Fig. [1] shows a schematic representation of several sim- 
ple hydrogen vacancy configurations on Pd(lll) consid- 
ered in the present work. A double vacancy can be con- 
sidered as a single H atom (red circles in Fig. [1]) mov- 
ing in the triangular area formed by 3Vh. As noted by 
Mitsui et al, this H atom is much more mobile than all 
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FIG. 1: (color online) Hydrogen vacancy configurations: (a) 
fully hydrogenated Pd(lll) surface; (b) single vacancy, IV; 
(c) classical model for a divacancy, 2Vc; (d) ground-state 
configuration of the quantum mechanical divacancy, 2Vq; (e) 
trivacancy centered around a hep site, (3Vh)- Pd atoms are 
represented by large gold circles and H atoms by smaller black 
circles. The white circles represent empty fee sites while the 
red circle represents a H atom that moves in the triangu- 
lar area formed by 3Vh. (f) An STM image of a divacancy. 
Reprinted from Ref. [3 with permission. 



other H atoms that form 3Vh.^^'^^ Previous models that 
treated H atoms classically considered only 2Vc where 
the (red) H atom hops among three fee sites via thermal 
diffusion. As we will show later, however, full quantum 
treatment of the mobile H atom reveals that the ground- 
state configuration for divacancy is 2Vq where the central 
H atom occupies the hep site rather than fee sites. Ex- 
change with the H atoms forming the walls of the triangle 
of 3Vh requires hopping near top Pd sites and/or close 
to other H atoms, a process with a higher barrier and 
therefore much slower than the motion of the H atom in 
the vacancy. Hence, the problem of determining the min- 
imum energy configuration of a divacancy is equivalent 
to solving the Schrodinger equation for a single hydrogen 
atom moving in the external potential of hydrogen and 
palladium atoms having the 3Vh configuration. 

Within the Born-Oppenheimer approximation^^, the 
Schrodinger equation for the motion of a hydrogen atom 
with mass M is 

- ^^'^(^) + U{rMr) = E^{v). (1) 



The 3-dimensional potential energy surface (PES), C/(r), 
contains the contributions to the total energy from all Pd 
atoms in the slab, all hydrogen atoms adsorbed on the 
surface, and the hydrogen atom in question at r. When 
only the position of this hydrogen atom is varied, U{r) 
becomes, in effect, the potential energy surface for its 
motion. 

The PES was mapped out by calculating the adsorp- 
tion energy of the hydrogen atom placed at different po- 
sitions over a 3Vh defect. We formed the defect within 
a periodic 3x3 surface unit cell, which is large enough 
to prevent significant interaction with periodic images in 
neighboring cells. The potential energy surface was sam- 
pled according to the importance of each region: sample 
points were densely distributed near points of interest 
such as the fee sites, the hep sites, and along the path- 
way between them, while points far from these locations 
were sampled less densely. Altogether over 3500 energy 
points were evaluated, distributed within about 60 planes 
parallel to the surface and separated by 0.1 A. 

All ab initio total-energy calculations and geome- 
try optimizations are performed within density func- 
tional theory (DFT)^^ using ultrasoft pseudopoten- 
tials (USPP)^^ as implemented by Kresse et. alM^ 
Exchange-correlation effects were treated within the 
local-density approximation (LDA) ^^i^^ . The electron 
wave functions were expanded in a plane-wave basis set 
with cutoff energy 250 eV. The resulting equilibrium lat- 
tice constant of bulk Pd was 3.86 A in excellent agree- 
ment with the experimental value of 3.89 A^^. We used a 
standard supercell technique to model the Pd(lll) sur- 
face as slabs of four Pd layers separated by 12 A of 
vacuum. The surface Brillouin zone was sampled us- 
ing a 6x6 Monkhorst-Pack /c-point set. Structural op- 
timization was performed until the energy difference be- 
tween successive steps was less than 1 meV. After ap- 
plying a well-known correction for the LDA binding en- 
ergy of H2r^, we obtained the hydrogen adsorption energy 
^ad=0.46 eV, in excellent agreement with the experimen- 
tal value of 0.45 eV— . Based on tests performed with 
thicker slabs and larger vacuum gaps, we estimate the 
uncertainty in this energy to be roughly 0.02 eV. Tests 
using the generalized-gradient approximation (GGA)^^^ 
show that £^ad values for different sites shift rigidly. This 
would have no effect on our main results below, because 
the PES depends only on differences in adsorption ener- 
gies, which are unchanged. 

We solved the Schrodinger equation ^ numerically in 
momentum space. Although we are modeling a hydro- 
gen atom moving in an isolated 3Vh defect, we nonethe- 
less use periodic boundary conditions in order to take 
advantage of Bloch's theorem and the convenience of us- 
ing plane waves as a basis. To minimize computational 
requirements, we used a unit cell with the full Csy sym- 
metry of the underlying lattice. 

The wave functions for the hydrogen atom are ex- 
panded in a plane-wave basis with a cutoff energy of 
0.3 eV, which is equivalent to a cutoff of roughly 550 eV 
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FIG. 2: (color online) Minimum potential energy surface for 
a hydrogen atom in a divacancy. The minimum potential 
energy in the z-direction for all points in the xy-plane of a 
single classical H atom moving in a three site cluster. Large 
and small white circles representing top layer Pd atoms and 
neighboring H atoms, respectively, are drawn as visual guides. 
The color bar shows corresponding potential values in eV. 

for electronic band-structure calculations. Our tests 
show that the eigenvalues of the hydrogen- atom wave 
functions are well converged with this cutoff. Since we 
simulate an isolated 3Vh defect with a large unit cell, the 
Brillouin zone is very small, and hence k point sampling 
is not important. The eigenvalues for different k points 
exhibit very little dispersion, and therefore we use just 
a single /c-point at the zone center. Finally, the wave 
functions are well confined inside the potential well, with 
no significant tailing outside the boundaries of the po- 
tential well, indicating that the 3x3 surface unit cell is 
adequate. 



III. RESULTS 

Fig. [2] shows the minimum PES for a hydrogen atom in 
a 3Vh defect. Two important conclusions can be drawn 
from this energy surface. First, the only low-energy path- 
ways available to the hydrogen atom are those connecting 
the three fee sites to central hep site. Second, there are 
no low-energy pathways available to allow the surround- 
ing hydrogen atoms to move into the vacancy region. The 
region of low potential (blue and green) is surrounded by 
neighboring H atoms and by construction possesses Csy 
symmetry about the central hep site. The potential floor 
(blue) region has three branches extending from the cen- 
tral hep site to fee sites, where the potential reaches its 
minimum. The local minimum at the hep site is 20 meV 
higher than at the fee sites. 

The results of our calculations for the low-energy eigen- 
states are summarized in Table IH and their wavefunctions 
are shown in Fig. [3l Of particular interest is the ground 
state (1^1 ), which is localized at the hep site. We note 
that this is not the minimum of the potential-energy sur- 
face (which occurs at the fee sites). If the H atom had 
been treated classically, the ground state configuration 



TABLE I: The energies of eigenstates for a H atom in a diva- 
cancy. The potential energies are measured from the bottom 
of the potential. Difference in energies are measured from 
those of the ground state (lAi). States 2Ai, lEi, and IE2 
are considered to be degenerate. All energy values are given 
in meV. 



Index 


state 


Etot 


^kin 


Epot 


A^tot 


A^kin 


A^pot 


1 


lAi 


177 


68 


109 








2 


2Ai 


191 


72 


119 


14 


4 


10 


3 


lEi 


192 


71 


121 


15 


3 


12 


4 


IE2 


193 


71 


122 


16 


3 


13 


5 


2Ei 


245 


78 


167 


68 


10 


58 




FIG. 3: (color online) Wave functions of a H atom in a di- 
vacancy. Isosurface of the real component of the wave func- 
tions of a hydrogen atom in a structure used to simulate a 
2Vq vacancy defect that contains three fee sites and one hep 
site. The isovalues on red surfaces are positive while the val- 
ues on blue surfaces are negative. All imaginary components 
are zero, (a) The ground state, The red surface con- 

tains 99.3% of the total integrated probability (TIP), (b) The 
first excited state, 02^1. Three red surfaces contain 99.1% 
of the TIP. (c) One of the first doublet excited states, (I)iei • 
The red surface contains 58.7% of the TIP while two blue 
surfaces contain 40.3%. (d) The other first doublet excited 
states, (l)iE2- Each surface contains 49.6% of the TIP. (e) One 
of the second doublet excited states, (I)2Ei- The red surface 
contains 61.4% of the TIP while two blue surfaces contain 
37.5%. (f) One of the mixed states delocalized over an fee 
site: ^fcc = ^(V2(/)2Ai - (piEi - V^(t)iE2)- The red surface 
contains 99.4% of the TIP. 
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would have the H atom adsorbed at the fee site, corre- 
sponding to the defect configuration 2Vc. 

To understand better the origin of this quantum effect 
we hst separately, in Table U the kinetic- and potential- 
energy contributions to the total energy. Our result in- 
dicates that the quantum-mechanical shift is driven by 
both the kinetic and potential energy of the hydrogen- 
atom wave function. Even though the minimum poten- 
tial at the hep site is slightly higher than that of three 
fee sites, the potential well at the hep site is wider than 
those at the fee sites. By localizing around the cen- 
tral hep site, the ground-state wave function can suffi- 
ciently lower its potential energy (by reducing the pen- 
etration of wave function into the region of higher po- 
tential) and kinetic energy (by reducing the gradient of 
the wave function) to compensate the small increase of 
potential around the hep site. Hence, the fully quan- 
tum mechanical ground-state geometry for the hydro- 
gen divacancy is not the classical configuration 2Vc of 
Fig. 1(c), but rather the quantum configuration 2Vq il- 
lustrated in Fig. 1(d). Indeed, in order to realize a con- 
figuration similar to the classical 2Vc, the quantum wave 
function must be localized on one of the fee sites as 
shown Fig. [3](e). This requires constructing, at the cost 
of 14 meV, a linear combination three degenerate excited 
states: ^fcc = 4^(V2^2Ai - (l^iEi - V30i£;2)- 



IV. DISCUSSION 

One prediction that arises from this fully quantum- 
mechanical calculation is that the ground-state of the 
hydrogen quantum divacancy 2Vq will have intrinsically 
triangular symmetry, suggesting that STM images should 
appear triangular as well. Because it is locally insulating, 
an adsorbed hydrogen atom screens the tunneling current 
between the STM tip and the metal substrate. Thus, a 
site occupied by an hydrogen atom has a low apparent 
height (dark in a standard STM representation) while 
a vacancy appears high (bright) As a result, a diva- 
cancy is expected to appear as a three-lobed object with 
triangular symmetry. Recent STM images, reproduced 



in Fig. mf), confirm this behavior— li^. Alternative ex- 
planations based on classical thermal diffusion^^^^^, al- 
though consistent with the observed symmetry, are not 
necessary to explain such images. For the trivacancy and 
larger multi-vacancy defects, the quantum effects mani- 
fested in the divacancy are suppressed and the delocalized 
H atoms behave more like classical particles. In these 
cases, thermal diffusion provides a reasonable explana- 
tion for the observed STM imagery^^'^^. 

Our results also suggest that future studies involving 
divacancies of H atoms on metal (111) surfaces need to in- 
clude the quantum configuration 2Vq shown in Fig. [IJd) 
as one of the viable stable configurations even if the total 
energy of the structure when H atoms are treated classi- 
cally is somewhat higher. 

V. CONCLUSIONS 

We have studied the quantum nature of hydrogen 
atoms on Pd(lll) surface by solving the Schrodinger 
equation for an H atom moving in static potential en- 
ergy surface determined from first-principles density- 
functional theory calculations. We find that a H atom 
in a divacancy defect experiences significant quantum ef- 
fects, with the result that its wave functions are extended 
over large portion of the vacancy. The ground-state H 
wave function is centered at an hep site rather than the 
fee site occupied by classical H atoms. As a result, the di- 
vacancy should exhibit a triangular geometry with three- 
fold symmetry, consistent with recent experiments. 
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